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The static dielectric properties of dense ionic fluids, e.g., room temperature ionic liquids (RTILs) 
and inorganic fused salts, are investigated on different length scales by means of grandcanonical 
Monte Carlo simulations. A generally applicable scheme is developed which allows one to ap¬ 
proximately decompose the electric susceptibility of dense ionic fluids into the orientation and the 
distortion polarization contribution. It is shown that at long range the well-known plasma-like per¬ 
fect screening behavior occurs, which corresponds to a diverging distortion susceptibility, whereas 
at short range orientation polarization dominates, which coincides with that of a dipolar fluid of 
attached cation-anion pairs. This observation suggests that the recently debated interpretation of 
RTILs as dilute electrolyte solutions might not be simply a yes-no-question but it might depend on 
the considered length scale. 


I. INTRODUCTION 

Room-temperature ionic liquids (RTILs) have become 
an extensively investigated topic in recent years. RTILs 
are salts with melting points at about room temperature, 
in contrast to inorganic salts such as NaCl with a melting 
point of 1081 K [l| . The first report of an RTIL is dated to 
1914, when the compound ethylammonium nitrate with a 
melting temperature of 285 K was synthesized [ 2 }. Since 
the 1990s interest in RTILs has grown greatly due to 
the discovery of a huge number of air- and water-stable 
compounds (see, e.g., Ref. Q and the references therein). 

Besides the strikingly low melting point, RTILs ex¬ 
hibit a number of remarkable general, i.e., material- 
independent, physical properties such as an extremely 
low vapour pressure Hi and a high viscosity at ambi¬ 
ent temperature — the lowest viscosity of RTILs at 298 K 
observed by now, 21 cP for [C 2 mim][N(CN) 2 ] Hj, is more 
than twenty times that of water. All these properties 
are related to the combination of strong electrostatic and 
highly anisotropic steric forces. 

An interesting feature of the RTILs for applications 
is the high sensitivity of some physical properties on 
the choice of the cation and the anion. For example 
the viscosity of [C 4 mim][BF 4 ] is only half of that of 
[C 4 mim][PFg] even though these RTILs differ only in the 
size of the (spherical) anions [T]. Therefore exploration of 
RTILs on the molecular level is crucial, in order to being 
able to “design” salts with prescribed properties. Com¬ 
puter simulations are frequently used as a tool, but the 
challenge is to find an adequate description of the inter¬ 
action potentials. Such interaction potentials are com¬ 
plicated and depend upon many parameters (e.g., bond, 
angle and torsion parameters as well as van der Waals 
terms Q). However, the aim of the present work is not 
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to model a specific RTIL with all its peculiarities but to 
investigate ionic fluids, i.e., RTILs and fused salts, from a 
general point of view which considers them mostly as flu¬ 
ids being composed of ions at high densities. Hence, since 
general features are expected to occur in all systems, one 
can focus on the simple system of equally-sized charged 
hard spheres, which is commonly called the restricted 
primitive model (RPM). 

As RTILs are used as solvents in chemical studies, their 
polarity is one of the most important characteristics be¬ 
cause it describes the global solvation capability of the 
solvent. In order to interpret surface force apparatus 
(SFA) measurements, it has been proposed recently Q 
to view pure ionic liquids as dilute electrolyte solutions 
with a few mobile ions in an effective solvent made of 
temporarily paired ions. Whereas these SFA data have 
been doubted mm, the interpretation of RTILs as 
dilute electrolyte solutions has attracted some interest 
0 ]. For non-conducting fluids the static dielectric con¬ 
stant e is well defined and precisely measurable (see, 
e.g., Ref. 0 , 111 ). However, conducting fluids, such as 
ionic fluids, are well-known to perfectly screen external 
charges at long ranges 0,0 so that the correspond¬ 
ing static dielectric constant is infinitely large. On the 
other hand, it is also well-known that cations and an¬ 
ions form a characteristic alternating pair structure at 
short ranges fl5l fl6| , which, when considering neighbor¬ 
ing cation-anion pairs, is analogous to the charge separa¬ 
tion inside dipolar molecules. Therefore, it appears that 
the dielectric properties of dense ionic fluids depend on 
the length scale. 

In the present work the analysis is based on the re¬ 
sponse of bulk ionic fluids to a static nonuniform electric 
field which spatially varies on a particular length scale. 
It has to be stressed that here the focus is on nonuni¬ 
form static electric fields, the response onto which is de¬ 
scribed in terms of the wavenumber-dependent dielectric 
function e(k), whereas numerous experimental studies 
consider uniform time-dependent fields, which give rise 
to the frequency-dependent dielectric function e(w) (see, 
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e.g., [T7M20|l Both quantities, e(k) and e(uj), are not eas¬ 
ily related because the former describes the equilibrium 
structure whereas the latter quantifies the dynamics in 
uniform electric fields. 

The dielectric properties of a substance can be in¬ 
terpreted in terms of two well-known mechanisms a 
Orientation polarization refers to the rotation of dipolar 
moments upon keeping the magnitude constant, whereas 
distortion polarization describes the change of magnitude 
of dipolar moments at constant orientations. The accord¬ 
ing orientation susceptibility Xori(fc) and distortion sus¬ 
ceptibility Xdis (k) both contribute to the total electric 
susceptibility x(fc) = e(k) - 1 = Xori (k) + Xdis(k). How¬ 
ever, in order to infer the dominant polarization mech¬ 
anism as a function of the wave number k, one has to 
somehow determine the decomposition of the observable 
x(k) = e{k) — 1 into the orientation and the distortion 
contribution. Perfect screening in ionic fluids and the ac¬ 
cording divergence x(.k) —> 00 in the limit k^O [HIIII 
corresponds to a dominating distortion polarization in 
the long-wavelength limit. This long-range behavior of 
ionic fluids is in sharp contrast to that of dipolar fluids, 
whose electric susceptibility y (k) attains a finite limit 
y(0) as k -» 0. 

In order to achieve the decomposition of the electric 
susceptibility x(k) of an ionic fluid into the orientation 
and the distortion susceptibility Xori (k) and Xdis(fc), re¬ 
spectively, the following approach is proposed: In addi¬ 
tion to determine the electric susceptibility x(^) of the 
ionic fluid composed of cations and anions, a correspond¬ 
ing dipolar fluid is considered whose particles are over¬ 
all charge-neutral dumbbells formed by gluing together 
pairs of cations and anions of the ionic fluid. This dipo¬ 
lar fluid does not exhibit distortion polarization, i.e., 
Xdis(fc) = 0, but only pure orientation polarization, i.e., 
x(k) = Xori (k), and the latter can be expected to be sim¬ 
ilar to Xori (k) of the corresponding ionic fluid at short 
ranges. The applicability of this approach is not re¬ 
stricted to the RPM investigated here but it can be used 
for any ionic fluid model. 

A brief description of the theoretical methods as well 
as the considered models is given in Sec. HU The results of 
actual calculations are discussed in Sec. m from which 
general conclusions on the dielectric properties of ionic 
fluids are drawn in Sec. m 


II. MODELS AND METHODS 
A. Restricted primitive and dumbbell model 

The model used to represent the ionic fluid is the three- 
dimensional restricted primitive model (RPM), i.e., a 
collection of N/2 positively and N/2 negatively charged 
hard spheres of equal diameter a and equal absolute va¬ 
lencies \z+\ = \z-\ (see Fig. [T]) a)). The interaction poten¬ 
tial between two ions of species i and j with i, j € {+, —} 



FIG. 1. Sketch of (a) the restricted primitive model (RPM) 
and (b) the corresponding dumbbell model (DM) with parti¬ 
cles being constructed by pairwise gluing together oppositely 
charged hard spheres of the RPM. 

at positions r, and tv,-, respectively, can be described as 

{ ZiZjl B , , . 

I'D - rj\ ' Tl Vj (1) 

00 , |n - rj\ < a 

with the vacuum Bjerrum length Ib = /3e 2 /(47T£o), where 
e is the elemenatry charge, (3 = 1/(&bT) denotes the 
inverse temperature, and £0 is the vacuum permittivity. 

In this work, the dielectric properties of the RPM (see 
Fig. [Ua)) are compared with those of a corresponding 
dipolar fluid model whose particles are composed of one 
cation and one anion of the RPM glued together (see 
Fig. m)- The particles within this dumbbell model 
(DM) possess three positional and two orientational de¬ 
grees of freedom, in contrast to six translational degrees 
of freedom of a pair of ions within the RPM. Obviously, 
all configurations of N/2 dumbbells correspond to possi¬ 
ble configurations of the RPM with N/2 positive and N/2 
negative hard spheres, but not all configurations of the 
RPM can be realized within the DM. The interaction en¬ 
ergy between two dumbbell particles is given by the sum 
of contributions Eq. © of the constituent charged hard 
spheres. 

B. Dielectric properties 

In this work the dielectric properties of ionic and dipo¬ 
lar fluids are studied by analyzing the linear response 
of the models introduced in the previous Subsec. Ill Al in 
the presence of a weak static nonuniform external elec¬ 
tric field. Within the linear response regime the dielec¬ 
tric properties are given by the dielectric function tensor 
^ (k) or, equivalently, by the electric susceptibility ten¬ 
sor *x/{k) = ^{k) — 1. However, since electrostatic fields 
are purely longitudinal due to Faraday’s law, the polar¬ 
ization field is also purely longitudinal in isotropic fluids. 
Therefore, only the longitudinal components (parallel to 
the wave vector k ) 

,(*):= = and 

*(*):= X|| (*) = AAAiA =e( fc)-i (2) 
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are of relevance here. 

It can be shown [22j that the longitudinal dielectric 
function e(k) is related to the charge-charge structure 
factor 


s zz (k) = — {^2 z m exp(-ik • r m )^ z n exp(ik ■ r n ) 


(3) 


via 


e{k) 




(4) 


In the present work the averaging in Eq. m is obtained 
by means of grandcanonical Monte Carlo simulations of 
the RPM and the DM. 


C. Grandcanonical Monte Carlo simulations 

In order to determine the charge-charge structure fac¬ 
tor S zz (k) in Eq. ([3]), grandcanonical Monte Carlo simu¬ 
lations with Metropolis sampling [23l - l25| were performed 
in a cubic box of side length E 1 / 3 = 10er with periodic 
boundary conditions applying Ewald’s method [241426} . 
In a first step, quick simulation runs were used to de¬ 
termine the relation between the chemical potential and 
the mean density N/V for fixed temperature. For calcu¬ 
lation of the charge-charge structure factor S zz {k) typi¬ 
cally 10 6 — 10 7 relaxation and 10' — 10 8 evaluation steps 
were used. A Monte Carlo step involved a random selec¬ 
tion of one of the following possible moves: 

• translation of an ion (RPM) or of a dumbbell (DM) 
inside a small cube-like environment of size d 3 with 
probability P tr ans, 

• insertion of a cation and an anion (RPM) or of a 
dumbbell (DM) with probability (1 — Ptrans)/2, or 

• removal of a cation and an anion (RPM) or of a 
dumbbell (DM) with probability (1 — Ptrans)/2. 

Insertions and removals of charge-neutral entities guaran¬ 
tee global charge neutrality during the whole simulation. 
Values of Ttrans and d were chosen in the way that the re¬ 
sulting rate of acceptance of the trial states did not drop 
below 30%. The Monte Carlo code has been validated 
by comparison to results obtained using the molecular 
dynamics package ESPResSo Sill. 


III. RESULTS AND DISCUSSION 

In the present work the RPM (Sec. Ill Al) is consid¬ 
ered as a representative ionic fluid. Its dielectric func¬ 
tion e(k) is determined by using Eq. Q via calculat¬ 
ing the charge-charge structure factor S zz (k), Eq. J3|, 
by means of grandcanonical Monte Carlo simulations 
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FIG. 2. Charge-charge structure factor of (a) the RPM 
(Sec. Ill Al) and (b) the DM (Sec. Ill Al) for the tempera¬ 
ture T* = 1 and various values of the packing fraction 
77 6 {0.08,0.14, 0.23,0.3}. 


(see Subsec. Hi Cl) . In the following discussion thermody¬ 
namic states are considered with the packing fractions 
rj = irNcr 3 /(6V) £ [0.08,0.33] and the temperatures 
T* = <j/Ib £ [1/3,1]. These conditions correspond to the 
plasma parameters T = 2?7 1 / 3 /T* £ [0.86,4.15]. Since 
the temperatures T* are well above the critical tempera¬ 
ture T* < 0.1 of the RPM [2l|, no influence of the vapor- 
liquid phase transition is expected to occur. Moreover, 
T* is also well above the temperature range T* < 1/4, 
where significant Bjerrum ion-pairing is expected to oc¬ 
cur [29[ , which renders the RPM a dipolar fluid already in 
the gas phase. The high temperatures T* used here allow 
one to study the crossover from plasma-like to dipolar- 
fluid-like behavior. The general trend of that crossover 
for varying 77 and T* is captured by a simple approxima¬ 
tive calculation described later. 

Figures a) and (b) display the charge-charge struc¬ 
ture factors S zz (k) of the RPM ionic fluid ISec. lII Al) and 
of the DM ISec. Ill Al . respectively. The calculation was 
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performed several (5-10) times for each packing fraction 
separately, both for the RPM and the DM. This proce¬ 
dure allows one to calculate the mean value of S zz (k ) 
as well as its standard deviation shown as error bars in 
Figs. [2{a) and (b). It is apparent that S zz (k) of the 
ionic fluid is sensitive to density changes only for very 
low packing fractions 77 . The slightly more pronounced 
oscillations of S zz (k) within the DM (Fig. [^b)) as com¬ 
pared to those within the RPM (Fig. [2](a)) are perhaps 
an artifact of the DM, within which two charged hard 
spheres are kept exactly at a distance er, whereas the 
principal peak of the cation-anion pair distribution func¬ 
tion within the RPM has a finite width [3Cj. 

Using Eq. Fig. [3] displays the inverse dielectric 
functions 1 /e(fc) of the RPM and of the DM. The perfect 
screening condition |lj|[ 22 j implies l/e(fe) —> 0 for fc —> 0 
for an ionic fluid due to the asymptotic behavior of the 
charge-charge structure factor | 22 | 


S zz (k) 


k ^>0 k 2 
~ 'k 2 


(5) 


with the Debye length 1 /k being given by k 2 = AttIbN/V, 
which is shown in the insets of Figs. Efc) and[3[d). Unfor¬ 
tunately, the error bars of the three leftmost data points 
for 77 = 0.23 (Fig. [31(c)) and of the four leftmost data 
points for 77 = 0.3 (Fig. (3[d)) in the plots of the inverse 
dielectric function l/e(k) of the RPM are too large to 
conclusively infer the functional form in the long range 
limit k —>• 0 at these packing fractions. However the 
error bars of the RPM data in the interesting range of 
medium and short separations, |fc|er > 27 r, and of the 
DM data are small. Since the DM can exhibit only ori¬ 
entation polarization, its dielectric function e(k —> 0 ) ap¬ 
proaches a finite value of the “dielectric constant” e( 0 ), 
which increases upon increasing the packing fraction 77 
(see Fig. [3]). 

In order to determine the dominant polarization mech¬ 
anism in an ionic fluid the inverse dielectric functions 
1 /e(fc) of the RPM and of the DM in Fig. [3] are con¬ 
verted into the electric susceptibilities y(fc) = e(k) — 1 . 
Figure U shows that the electric susceptibility y RPM (fc) = 
X RPM (fc) + X RPM (fc) of the RPM, which comprises a con¬ 
tribution X' R ; (fc) due to orientation polarization and 
a contribution X^ PM (fe) due to distortion polarization, 
almost coincides in the range |fc|er > 27r with the elec¬ 
tric susceptibility x DM (fc) = Xor^(^) °f the DM, which, 
by construction, exhibits only orientation polarization 
with susceptibility xS^^)- Since for sufficiently large 
packing fractions 77 it can be expected that the orien¬ 
tation susceptibility Xori M (^) °f the RPM is identical 
to the electric susceptibility Xor^W °f the DM, i.e., 
XotT M (^) = Xot^(^) f° r wavenumbers fc, one can in¬ 
fer the orientation and the distortion susceptibility of the 
RPM separately: 


RPM 

Xori 


(fc)=x DM (fc), 
x£™(fc) = x RPM (fc)-x DM (fc). 


( 6 ) 



FIG. 3. Inverse dielectric functions l/e(fe) of the RPM ionic 
fluid (red solid lines, see Sec. Ill Al) and of the DM dipolar 
fluid (blue dashed lines, see Sec. Ill Al at temperature T* = 1 
and packing fractions r/ € {0.08,0.14,0.23,0.3}. Whereas 
l/e(fc —>• 0) becomes small (i.e., e{k —»• 0) becomes large) 
for the ionic fluid, the dielectric function of the dipolar fluid 
approaches a finite value e(0) for small wavenumbers. The in¬ 
set in panels (c) and (d) compares the charge-charge structure 
factor S zz (k) obtained by means of Monte Carlo simulations 
(red circles) with the asymptotic behavior Eq. © (black solid 
line). 


Figure [5] clearly indicates that, at T* = 1, orien¬ 
tation polarization is the dominant mechanism of the 
RPM at sufficiently large wave numbers fc, whereas dis¬ 
tortion polarization is dominating at sufficiently small 
wave numbers fc. The weak distortion polarization, i.e., 
X RPM (fc) ~ 0 at large wave numbers fc (see Fig. [5]) can be 
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FIG. 4. Electric susceptibility x{k) = £ {k) — 1 of the RPM 
ionic fluid (red solid lines, see Sec. Ill All and of the DM dipolar 
fluid (blue dashed lines, see Sec. Ill All at temperature T* = 1 
and packing fractions r/ € {0.08,0.14,0.23,0.3}. Whereas 
x(k —> 0 ) diverges for the ionic fluid due to the perfect screen¬ 
ing of the external charges, the susceptibility of the dipolar 
fluid approaches a finite value x( 0 ) for small wavenumbers. 
At large wavenumbers \k\a > 2n the electric susceptibility 
x{k) = Xori(fc) + Xdis(fc) of the RPM almost coincides with 
that of the DM, which, by construction, possesses only ori¬ 
entation polarization whose electric susceptibility can be ex¬ 
pected to be close to Xori(fc) of the RPM (see Sec. Ill At . 


attributed to the impenetrable hard cores of the ions of 
the RPM. The crossover wavenumber between distortion- 
dominated and orientation-dominated polarization in¬ 
creases with packing fraction 77. There is a certain inter¬ 
val of wavenumbers |fe|cr « 4 ... 5 with XdJ M (^) < 0 ( see 



FIG. 5. Distortion susceptibilities Xdis M (k) (solid lines) 
and orientation susceptibilities Xari(k) (dashed lines) of 
the RPM at temperature T* = 1 for packing fractions rj £ 
{0.08,0.14,0.23,0.3}. Perfect screening corresponds to the 
divergence of XdJ M (&) i n the long wavelength limit |fc| —► 0 . 
The region of negative distortion susceptibility, Xdis M (k) < 0, 
can be interpreted as overscreening. At sufficiently large wave 
numbers orientation polarization dominates over distortion 
polarization with the crossover wave numbers increasing with 
the packing fraction rj (see Fig. ED). 


Fig. [5j), which indicates overscreening, which is caused 
by steric effects, i.e., by the hard ion cores, too. The 
distortion susceptibility of the RPM diverges for k —> 0 
according to Xdi™(fc) ~ 1 /fc 2 , which corresponds to the 
perfect screening property of plasmas [laQ- Hence, it 
has been shown here, that the RPM ionic fluid exhibits 
dielectric properties similar to a dipolar fluid at short 
range, whereas it behaves plasma-like at long range. 

Figure [ 6 ] displays the crossover wave numbers 
k x (r),T*)<j (dots •), defined by equal orientation and 
distortion susceptibilities, Xod M (^ X ) = Xdis M (^ X )i as 
functions of the packing fraction 77 and of the temper¬ 
ature T*. The RPM exhibits plasma-like behavior for 
k < k x (r/,T*) and dipolar-fluid-like behavior for k > 
k x ( 77 , T*). The general trend is that of increasing values 
of k x (rj,T*)cr upon increasing the packing fraction 77 or 
decreasing the temperature T*. However, it can be ex¬ 
pected that k x (r],T*) < 2-K/a for any set of parameters 
( 77 , T*) since an external electric field oscillating with a 
wave length equal to the ion diameter a (see the thin ver¬ 
tical dashed line in Fig. [G]) cannot lead to distortion po- 
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FIG. 6. Crossover wave numbers k x (j],T*)a of the RPM 
between plasma-like behavior for k < k x (r/,T*) and dipolar- 
fluid-like behavior for k > k x (r/,T*) as functions of the pack¬ 
ing fraction 77 and of the temperature T*. The crossover wave 
number k x (r/,T*) is defined by equal orientation and distor¬ 
tion susceptibilities, XotT M (& X ) = Xdis M (^ X ) ( see Eq. ©). 
The dots (•) correspond to simulation data whereas the 
dashed lines are obtained by means of the simple approxi¬ 
mative calculation described in the main text. The latter 
captures the correct trends of increasing values fc x ( 77 , T*)<j 
upon increasing the packing fraction 77 or decreasing the tem¬ 
perature T*. The thin vertical dashed line marks ka = 2n, 
which corresponds to an upper bound of k x (r/,T*)a. 

larization, i.e., the crossover X' R T M (^ X ) = Xdi™(^ X ) h &s 
to occur at some smaller wave number. The dashed lines 
in Fig. [6] correspond to an approximation of fc x (? 7 ,T*)a 
with x RPM (^) Eq. ([ 6 ]) being approximated within the 
mean spherical approximation (MSA) of the RPM [22| 
and with x DM (^) i n Eq. (0 being approximated by the 
expression of a single dumbbell particle in an external 
electric field. Whereas this simple approximation slightly 
overestimates the value of k x ( 77 , T*)a, the general trends 
of increasing values of k x (t),T*)<j upon increasing the 
packing fraction 77 or decreasing the temperature T* are 
captured correctly. Hence for dense ionic fluids, e.g., 
inorganic fused salts (77 ~ 0.5, T* « 1/30) or RTILs 
(77 « 0.5, T* ss 1/50), one can expect plasma-like be¬ 
havior in a very wide range of wave numbers k < k x 
with 27r/fc x corresponding almost to the size of the ions. 


IV. CONCLUSIONS 

The main observation of the previous Sec. mu that 
the RPM ionic fluids exhibits dielectric properties simi¬ 


lar to a dipolar fluid at short range whereas it behaves 
plasma-like at long range, does not hinge on any pecu¬ 
liar property of the RPM and can hence be expected 
to be made for other ionic fluids, too. Moreover, our ap¬ 
proach in Sec.EHto decompose the electric susceptibility 
x(fc) = Xori(fc) + Xdis(fc) into a contribution x'ori(fc) due 
to orientation polarization and a contribution Xdis{k) due 
to distortion polarization by introducing a correspond¬ 
ing dipolar fluid made of cation-anion compounds ap¬ 
plies to the case of other ionic fluids, too. In general, 
the structure of the cation-anion compounds comprising 
the dipolar fluid corresponding to an ionic fluid can be 
conjectured on the basis of, e.g., the pair distribution 
function, which is routinely calculated for numerous ionic 
fluid models (see, e.g., Ref. [El,HU). Hence, the method 
described in Sec. EH could be a useful tool to assess the 
bulk dielectric properties of general ionic fluid models. 


The results of Sec. EH are restricted to bulk ionic fluids 
and they cannot be applied quantitatively in the context 
of the discussion on the interpretation of SFA measure¬ 
ments in an RTIL environment , which is related 

to confined ionic fluids. The reason for this restriction 
is that the static dielectric function for non-uniform sys¬ 
tems is of the form e(fc, k') due to the absence of trans¬ 
lational symmetry. However, the qualitative picture oc¬ 
curring in Sec. Mil suggests, that strongly confined ionic 
fluids tend to behave as dipolar fluids whereas they pro¬ 
gressively exhibit plasma-like properties upon relaxing 
the confinement. This suggests that the recently de¬ 
bated interpretation of RTILs as dilute electrolyte so¬ 
lutions :9i, 113] might not be simply a yes-no-question but 
it might depend on the considered length scale. 


The finding in Sec. EH that the static dielectric prop¬ 
erties of ionic fluids depend on the length scale, may be 
considered complementary to the observation that the 
quantification of the polarity of RTILs dep ends on the 
intrinsic time scale of the measurement [ 20 ( 1 . 


In summary, based on grandcanonical Monte Carlo 
simulations of the restricted primitive model and a corre¬ 
sponding dumbbell model, it has been argued that bulk 
ionic fluids at small length scales are expected to behave 
as dipolar fluids, for which orientation polarization domi¬ 
nates, whereas at large length scales plasma-like behavior 
occurs, for which distortion polarization dominates. In 
the style of Ref. [2C|, one can conclude that the static 
dielectric properties of ionic fluids depend on the length 
scale on which they are looked at. 
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